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The dynamical critical behavior of a single directed line driven in a random medium near the 
depinning threshold is studied both analytically (by renormalization group) and numerically, in 
the context of a Flux Line in a Type-II superconductor with a bulk current J. In the absence of 
transverse fluctuations, the system reduces to recently studied models of interface depinning. In most 
cases, the presence of transverse fluctuations are found not to influence the critical exponents that 
describe longitudinal correlations. For a manifold with d = 4 — e internal dimensions, longitudinal 
fluctuations in an isotropic medium are described by a roughness exponent C,\\ = e/3 to all orders 
in e, and a dynamical exponent zy — 2 — 2e/9 + O(e^). Transverse fluctuations have a distinct and 
smaller roughness exponent (^x = C|| ~ d./2 for an isotropic medium. Furthermore, their relaxation 
is much slower, characterized by a dynamical exponent z± = zy + 1///, where v — 1/(2 — (^y) is the 
correlation length exponent. The predicted exponents agree well with numerical results for a flux 
line in three dimensions. As in the case of interface depinning models, anisotropy leads to additional 
universality classes. A nonzero Hall angle, which has no analogue in the interface models, also affects 
the critical behavior. 
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I. INTRODUCTION AND SUMMARY 

The study of dynamical critical phenomena associ- 
ated with the pinning-depinning transition in random 
media has become a subject of considerable interest in 
recent years. This is due to the importance of pinning 
in a wide variety of technologically important phenom- 
ena such as flux line (FL) motion in Type-II supercon- 
ductors, dynamics of interfaces (phase boundaries, inva- 
sion fronts, cracks, surface growth, to name a few), and 
charge-density wave (CDW) transport. These systems 
are characterized by a rough energy landscape due to the 
randomness in the medium. At zero temperature there 
are two distinct "phases" , distinguished by an order pa- 
rameter (henceforth called velocity) that measures the 
dynamic response, such as the average velocity for a FL, 
or current for a CDW. For small driving forces, the sys- 
tem is trapped by one of the many available metastable 
stationary states, and is "pinned" to the impurities in 
the medium. Critical behavior emerges as the station- 
ary states disappear, and the system starts moving with 
a nonzero velocity, when the driving force is increased 
above a threshold value. Extensive experimental [Q, the- 
oretical [^[jij, and simulation Q work has been done to 
understand the properties of this transition in CDW sys- 
tems. There are also numerous studies on the depinning 
of driven interfaces P-pl]|. A better theoretical under- 
standing of this dynamical phase transition was recently 



achieved, and critical exponents were calculated through 
an e-expansion for both CDW systems |^ and driven in- 
terfaces More recently, we performed similar cal- 
culations for the depinning of an elastic line in a bulk 
random medium, like a polymer in a gel network, a FL 
in a typc-II superconductor, or a screw dislocation in a 
crystal . In this article, we present a detailed report of 
our study on the dynamical critical behavior associated 
with the depinning of a FL, and in general on the de- 
pinning of directed manifolds in random media, through 
methods similar to those used for CDWs and interfaces. 

Specifically, let us consider the geometry of the FL 
shown in Fig. |^. The superconductor is subject to a 
magnetic field B = Bx along the x-axis, and a bulk su- 
percurrent J = Jz along the z-axis. A FL is oriented 
along B on the average, but deviates from a straight line 
due to impurities in the superconductor, which are rep- 
resented by a potential V{x,y,z). The conformations of 
the FL are described by R{x,t) = xx + r{x,t), where 
v{x,t) — y{x,t)y -\- z{x,t)z is a two component vector, 
lying in a plane normal to the magnetic field. The bulk 
current J drives the FL along the y-direction through the 
Lorentz Force = x B. (<I>o is the flux quantum.) 
If the bulk current is large enough, the FL drifts with 
an average velocity v. Due to the chiral nature of the 
supercurrents around the FL, v is in general not along 
the y-direction, but makes an angle cf) with the y-axis. 
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This is usually called the Hall angle, and although typi- 
cally small |13| , it can be significant near the depinning 
transition. 

It is more convenient to work with components of r 
that are parallel and perpendicular to v, i.e. 



r{x,t) = r||(a;,t)e|| + rj_{x,t)e±, 




(a) 




(b) 

FIG. 1. Geometry of the FL in a medium with impurities: 
(a) Three-dimensional geometry, (b) A cross section of the 
medium at fixed x. The average drift velocity v = ney malces 
an angle <j) with the j/— axis. 

where the unit vectors e|| and are rotated by d from 
the y- and z-axes respectively, as shown in Fig. |l|b. In 



Sec.0 we show that, under very general assumptions, 
the equation of motion for small deviations around a 
straight line, generalized to d-dimesional internal coor- 
dinates X <E 9?'', can be written as 

T^dtrn =Xii V^ri, + Ki2ylr±+F+fn (x, r(x, t)), (1.2a) 



(IT) ^dtr^^K2iS/lr. + K22ylr± + f±ix, r(x, t)), 



(1.2b) 



where rj is the viscosity the FL and F = ^qJ. The mod- 
uli Ka'j relate the elastic force to the local curvature and 
are in general nondiagonal for a sample with orientation- 
dependent core energy, or nonzero Hall angle (cf. Sec. ^). 
The random forces fa that arise from the impurity po- 
tential V are taken to have zero mean with correlations 

(/„(x,r)^(x',r'))-<5'^(x-x')A„^(r-r'), (1.3) 

where A is a function that decays rapidly for large values 
of its argument. (The indices a, 7, . . . = {|| , J-}.) 

Ignoring fluctuations of the FL transverse to the di- 
rection of average velocity, i.e. setting r± = 0, leads 
to an interface depinning model studied by Nattermann, 
Stepanow, Tang, and Leschhorn (NSTL) 0, and by 
Narayan and Fisher (NF ) Hence, the major dif- 

ference between Eqs.(L2) (henceforth called the "vector 
depinning model" ) and the previously studied "interface 
model" is the existence of transverse fluctuations, mak- 
ing the position of the line r a vector instead of a scalar 
"height" variable. The effects of such transverse fluc- 
tuations for large driving forces and average velocities, 
when the randomness in the medium can be approxi- 
mated as uncorrelated in space and time, were shown 
p^ , p^ to create a much richer dynamical phase diagram 
than the corresponding interface growth model, namely 
the Kardar-Parisi-Zhang (KPZ) equation ||l^. Then, the 
natural questions to ask are: How do these transverse 
fluctuations scale near the depinning threshold, and how 
do they influence the critical dynamics of longitudinal 
fluctuations? 

In order to make these questions more quantifiable, we 
consider the exponents that characterize the critical be- 
havior near the depinning transition. Let F(v) denote 
the driving force required to move the FL with a velocity 
V ~ we||. For small values of F = |F|, the line is pinned 
by the disorder in the medium. There is a threshold force 
Fc, such that the line moves with a nonzero average ve- 
locity V iS F > Fc For F slightly above threshold, 
we expect the average velocity to scale as 



AiF - F,) 
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(1.4) 



where /3 is the velocity exponent and A is a, nonuni- 
versal constant. Superposed on the steady advance of 
the line are rapid "jumps" as portions of the line de- 
pin from strong pinning centers. Such jumps are similar 
to avalanches in other slowly forced systems and have a 
power-law distribution in size, cut off at a characteristic 
correlation length ^. On approaching the threshold, ^ 
diverges as 
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(1.5) 



defining a correlation length exponent v. At length scales 
up to ^, the interface is self-affine, with correlations sat- 
isfying the dynamic scaling form 

([r||(x,t)-r||(O,0)]2) = |xp«lig||(t/|xril), (1.6) 

([r^(x,<)-r^(O,0)]2) = |x|2«-g^(V|xr-), (1.7) 

where C.a and Za are roughness and dynamic exponents, 
respectively. The scaling functions ga go to a constant as 
their arguments approach 0; C|| and C.± are the longitu- 
dinal and transverse wandering exponents of an instan- 
taneus line profile; and z± characterize scaling of re- 
laxation times of longitudinal and transverse modes with 
wave vector through Ta{q) ~ q~^°'. Beyond the length 
scale ^, regions move more or less independently of each 
other and the system is no longer critical. The behavior 
of the movin g li ne is described by the exponents calcu- 
lated earlier |l^,^ for time dependent noise. Ignoring 
any potential nonlinearities leads to a regular diffusion 
equation with white noise, for which the roughness and 
dynamic exponents are = C]^ = (2 — d)/2, z+ = 2. In 
the interface model, transverse fluctuations do not exist, 
thus, Cj_ and zj are not defined. 



Equations (1.2) can be analyzed using the formalism of 
Martin, Siggia, and Rose (MSR) . A renormalization 
group (RG) treatment of the "interface model", stud- 
ied by NSTL and NF §, indicates an upper critical 
dimension of dc = 4, and exponents in d = 4 — e dimen- 
sions, given to one-loop order as( = e/3 and z — 2— 2e/9. 
NSTL obtained this result by directly averaging the MSR 
generating functional Z, and calculating the renormaliza- 
tion of the force-force correlation function A(r). NF, on 
the other hand, expanded Z around a saddle point so- 
lution corresponding to a mean-field approximation ||ig| ] 
to Eqs.( |l.2| ) which involves temporal force- force correla- 
tions C{vt). They point out some of the deficiencies of 
conventional low- frequency analysis, and suggest that the 
roughness exponent is equal to e/3 to all orders in per- 
turbation theory. They also show that for two different 
classes of disordered systems, random-field and random- 
bond disorder, the zero temperature interface dynamics 
is essentially the same near threshold. Their argument 
remains valid for vector depinning, and our results will 
be applicable to both types of randomness. As we shall 
demonstrate in Section II] , the longitudinal exponents of 



the "vector" model are identical to those of the depinning 
interface, and given by 

C||=e/3, (1.8) 
0|l =2-2e/9 + 0(e2). (1.9) 

Other exponents are determined by exact exponent iden- 
tities from C|| and Z|| as 

1 3 



2-Cll 
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/3=(z||-C||)^^ = l-e/9 + 0(e2). 



(1.10) 
(1.11) 



Following the formalism of NF, we employ a pertur- 
bative expansion of the disorder-averaged MSR partition 
function around a mean-field solution for scalloped impu- 
rity potentials . We show that slightly above threshold, 
transverse fluctuations do not significantly affect the dy- 
namics of longitudinal fluctuations, apart from shifting 
the threshold force Fc. Specificall y, the ex ponents and 
exponent identities given in Eqs. (1.8-1.11) for d < dc 
are also correct for the vector depinning model. How- 
ever, transverse fluctuations turn out to scale differently, 
with 7^ (^11 and z± ^ zy. In particular, in an isotropic 
medium with Hall angle ^ = (Model A in Section ||), 
the renormalization of transverse temporal force-force 
correlations C±{vt) yields 



C-i 



5e 



(1.12) 



correct to all orders in e. The transverse dynamic expo- 
nent is given by an exact exponent identity: 



zj. = Z|| + -=4--+0(e2). 



(1.13) 



These conclusions can also be generalized to more than 
one transverse direction: the results do not depend on the 
number of transverse coordinates. For the FL (e — 3), 
the critical exponents arc then predicted to be 



C|| = 1, z\i w 4/3, V = 1, 
P ~ 1/3, a = 1/2, z± « 7/3. 



(1.14) 



This implies that in a type II superconductor driven 
slightly above threshold, flux lines are contained mostly 
in the plane normal to the current, up to the correlation 
length scale ^. This may have a noticeable effect on the 
dynamics of entanglement of flux lines near depinning. 
These results also rationalize the use of a "planar ap- 
proximation" in numerical simulations of FL depinning 



Another important consideration is the role of 
anisotropy in the bulk material. It was recently shown 
that anisotropy leads to new universality classes in inter- 
face depinning |^l|. We show that this happens as well 
for FL depinning, in an even richer fashion. The presence 
of a nonzero Hall angle affects the critical behavior in a 
manner similar to anisotr opy. These issues are discussed 
in more detail in Section VIH . 



The rest of the paper is organized as follows: In Sec- 
tion U, we derive the general form of the equation of mo- 
tion for a single FL, starting from a reparametrization 
inva rian t (RI) descpription of the FL dynamics. In Sec- 
III , we flrst establish the connection of Eqs. ( pT^ ) to 



tion 



the interface depinning problem for the simple case of an 
isotropic medium with zero Hall angle. We then examine 
the linear respo n se of the system to derive the exponent 
identities (p^) , (pTlll) , and dTT^ ), which are later shown 
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to be consistent with a formal RG treatment of the prob- 



lem in more general circumstances. In Section IV 



present the MSR formalism and expand the generating 
functional around a self-consistent saddle point solution, 
given by a mean-field theory. In Section we calcu- 
late response and connected correlation functions of the 
mean-field theory, which correspond to the bare propa- 
gators and vertex functions in a perturbative expansion. 
In Section VI, we determine critical exponents through 
an e-expansion near d = 4 dimensions, and in Sec. VII 



we compare these with numerical results obtained by di- 
rectly integrating the equations of motion. Finally, in 
Section VII] we discuss the physical significance of these 



results, the roles of nonlinear terms and anisotropy, and 
applicability of similar methods to related problems. 



II. EQUATIONS OF MOTION FOR A FL 

In this section we derive a phenomenological equation 
that describes the coarse-grained (in space and time) evo- 
lution of a single FL in a bulk type-II superconductor. 
The configuration of the FL at time t is described by 
i?(s, t), where s is an arbitrary parameter which we shall 
later equate to the a;— component of R. The equations 
of motion are obtained by balancing the "conservative" 
and "dynamical" forces. Conservative forces are derived 
from the energy functional and depend only on the in- 
stantaneous configuration R{s) of the FL. They include 
the elastic force, random forces due to the impurity po- 
tential V , and the Lorentz force due to the bulk current. 
Dynamical forces, on the other hand, depend explicitly 
on the local velocity of the FL and comprise the dissipa- 
tive and Magnus forces [p2| . 

For notational simplicity, we set the external magnetic 
field B along the x-axis and the the average velocity 
V along ey, suppressing the possible dependence of pa- 
rameters on the relative orientation of B and ey due to 
anisotropy in the underlying mate rial. Such complica- 
tions will be taken up later in Sec. VIII. An important 



consideration is the requirement that the equation of mo- 
tion be invariant under an arbitrary reparametrization 
R{s) R{s') of the curve. One such reparametriza- 
tion invariant quantity is the infinitesimal arclength dl = 
dsy/g, where g = dsR • dsR is the metric. The only phys- 
ically observable motions of the FL are orthogonal to the 
local unit tangent vector 



t 



-^dsR. 



Assuming that the FL motion is overdamped, the con- 
servative force Fx, which is derived below, is balanced 
by dynamical forces that are proportional to the local 
normal velocity = V ■ dtR = dtR — {dtR ■ i)t. (Here, 



titj projects any vector onto the local normal 



Vn- In general, there is an angle 4> (called the Hall an- 
gle) between the applied force and the velocity of the FL. 
Physically, this is due to the Magnus force which is or- 
thogonal to the velocity, and the Hall effect in the normal 
core of the FL ||2^ . The equation of motion can then be 
written as 



rjV ■ ^cos(t)dtR + sin(j){dtR) x i| Fj 



(2.1) 



To determine the conservative force Ft, consider the 
energy cost associated with a particular coarse-grained 
configuration R{s) of the FL in the absence of a bulk 
current, which is 



E 



i?(s) 



ds^/g- 



dsR ■ (T ■ dgR 

9 



V{R{s))\. (2.2) 



In the above equation, the symmetric tensor cr gives the 
core energy per unit length of the FL, and can be non- 
diagonal for an anisotropic sample. (Anharmonic contri- 
butions to the core energy can be ignored in a coarse- 
grained description and we will systematically keep only 
the leading order elastic terms.) The restoring force Fb 
is given by the energy cost of an infinitesimal virtual dis- 
placement SR{s). After some rearrangement, we arrive 
at 



dE 



j dsy^SR ■ V ■ |2cr • K - (< • (T • t) K 

+V{R)n~W^V{R)^ 
ds^5R-FB, (2.3) 



where k = g ■ d^R is the local curvature vector. To 
leading order, the random potential V{R) that multiplies 
K can be approximated by its spatial average, and elim- 
inated without loss of generality by choosing {V) = 0. 
/ = —WjiV{R) acts as a random force on each segment 
of the FL, whose correlations in general satisfy 



{UR)f^iR'))^A^^{R-R') 



(2.4) 



plane.) Dynamical forces are not necessarily parallel to 



For now, we do not restrict the form of A, apart from 
the reasonable expectation that it decays quickly beyond 
a characteristic impurity size a. When a bulk current 
J is present, the FL is also subject to a Lorentz force 
Fl = ^oJ X t, where $o is the flux quantum. Thus, the 
total conservative force acting on a section of the FL is 
given as 

Ft=v I 2 — _J + $0 J X i + / ^ . 

(2.5) 

For an isotropic sample in the extreme type-II limit, the 
Bardeen-Stephen model gives l23] 
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o'ij ~ dij 

77«<i>^/(2<V^n), 

where is the London penetration depth, is the co- 
herence length, and p„ , p^ are normal and Hall resistivi- 
ties of the non-superconducting core region, respectively. 
More general expressions for these phenomenological pa- 
rameters can be derived from a mesoscopic model based 
on a time-dependent Ginzburg-Landau theory pst . 

Equation ( |2.1[ ) is highly nonlinear and generalizes those 
of ref. 1^ to the three-dimensional and anisotropic case. 
We now pick {a;,e||,ex} as our coordinate axes, and x 



as the arbitrary parameter s, representing the FL as 
R{x,t) = XX + r||(a;,t)e|| -I- r±{x,t)e±. In this repre- 
sentation, g = 1 + {dxr\\)^ + {dxr±y, J — J||e|| + J±e±, 
/ = fxX + /||e|| + /_ie^, and 



V crx± 



<7x. 



After some rearrangement, and elimination of higher- 
order terms coming from the elastic energy of the FL, 
we obtain the following evolution equations for the com- 
ponents r|| and rj_: 



COS' 



-dtr\\ = [(2(T|| - ax) ~ 2ax tan(/)]9^r|| -I- [2ax ~ (2(Tj_ - (Tx) tan 0]9^r_L 



+^ {J±[l + {dxr\\)^] - J\\ [dxrudxrj^ - tan^V?]} 
v5 



tan(? 



■dxr\\dxri_ 



-I 



tan cp 



[I + (dx^)'] - fx 



tan (/ 



■dxrA 



(2.6a) 



cos I 



:dtr± = [2(Tx + (2cr|| - cr^) tan0]a^r|| -I- [{2aj_ - ax) + ay, tan0]9|r_L 



$0 

V9 



{J± [dxr\idxr± +ta.n(j)^] - Jy [l-f {dxr±)'^]} 



tan ( 



+f\l^[l + idxr^)^]+h 
v9 



1 tan (/) 

1 —Oxr\\axri_ 



V9 



fx 



« , tan(/> 
Oxr± H ^Oxrn 

V9 



(2.6b) 



These equations are clearly t oo c omplicated for an exhaustive analysis. However, it is possible to perform a gradient 



expasion of the RHS of Eqs.(|2.6|) when the fluctuations around the straight line are small, i.e. {dxTii)'^, {dxr±)'^ <C 1. 



In that case, Eqs.(p.6|) simplify to 



C0S( 



■dtvii = [{2aii — ax) — 2ax tan(/)]9fr[| + [2ax — {2a± — ax) ta,n(p]dxr± 
+^o{j± + J\i tan(/))+/||-/j_tan(/). 



(2.7a) 



COS I 



7dtr± = [2(Tx + {2a\\ - ax) tan0]9^r|| -I- [(2a j_ - ax) + a^ tan(/)]5^r_L 

+$o(^± tan(/) - J||)+/±-/|| tan(/), (2.7b) 



neglecting all terms of O ((S^rn)^, (dxr±)'^) or higher. 

So far, we have not enforced the condition that e|| points along the average velocity of the FL. This is satisfied by 
the self-consistency relation 

{dtri_) = 0. (2.8) 

In the small fluctuation limit where Eqs.( p/7| ) are valid, this condition is satisfied simply by setting J|| = J±ta,n(p. 
In order to study the scaling properties of this system in the framework of a field theory, we generalize the FL to a 
manifold with d-dimensional internal coordinates x € SR''. Further rearrangements, and addition of an infinitesimal 
external force £(x, t) in order to study response functions, lead to 

r]dt^ =ifiiV2r||+i^i2V2ri+F + /^|(x,r(x,t))+ei(x,t), (2.9a) 
r^dtri_ = K^iVln + K22Vlr^ + /l(x, r(x, t)) + e2(x, t), (2.9b) 
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where F ^ ^n. J?; + J"] = $o J, and 



K21 K22 J \sm(p coscj) J \ 2(Tx '2a± — (Tx 

( h \ = ( cos(/i -sin(?!)\ / /|| \ 
\f±J \sm(l) cos^ J \ J ' 

The correlations of the random forces satisfy 

(/„(x, r)^(x, r')) = ^''(x - x')A„^(r - r'). (2.9c) 



(Note that while both r and x are represented by bold 
characters, r remains two dimensional, while x has been 
promoted to a d-dimensional vector.) 

In the special case of an isotropic medium with = 0, 
the equations further reduce to 

jjStrii =ii'V2r|| +F + /||(x,r(x,i)) +£i(x,i), (2.10a) 

T]dtri_ = KVlr^ + /^(x, r(x, <)) + e2(x, <), (2.10b) 

where the correlations of the random forces satisfy 

(/„(x, r)^(x, r')) = S^y{x - x')A(|r - r'|). (2.10c) 

We shall henceforth refer to Eqs.( ^.ld ) as Model A. 
Anisotropy and/or a nonzero Hall angle changes the scal- 
ing properties of the critical region, and we sh all refer to 
this more general case, described by Eqs.(2^), as Model 
B. 



III. THE VECTOR DEPINNING MODEL 

In this section, we study some prop erties of the sys- 
tem described by Eqs.(|2.9|) and ( 2.10 ), in detail. Due 
to statistical translational symmetry in time t and in- 
ternal coordinates x, we use the real (x, t) and Fourier 
(q, uj) domains interchangeably when dealing with statis- 
tical averages. 

The vector depinning model differs from the CDW or 
interface problems due to the presence of transverse fluc- 
tuations r_L(x, i). It is sometimes useful to recast the 
equations such that r± appears as a function of ry rather 
than t. The asymmetry in r\\ and r± occurs because r|| 
almost always moves in the forward direction p^ , and 
therefore is a monotonous function of t. Thus, for any 
particular realization of the random force /(x, r), there 
is a unique point 7'^(x, ry) that is visited by the line 
for given coordinates (x, rii). The evolution of rj_ (x, r\\) 



A. Model A 



can be obtained schematically, by dividing Eq.(2.9b) by 
(Ua), as 



dr± 

dr\\ 



K2iVlr\\ + K22Vlri_ + fi 



(3.1) 



/fiiVjrii + K^2^lr^ + /|| + F 

We shall see that in most cases, the scaling properties 
of r_L in relati on t o can be obtained heuristically by 
inspecting Eq.(3.1). 



F irst of all, we establish the connection between 
Eq.( ^.ld ) and the interface depinning model for the spe- 
cial case of an isotropic system with (f) — Q (Model A). For 
a particular realization of randomness f (x, r), Eq.( p.lC| a) 
can be written as 



r^dtr\\ = ifV^rii + /' (x, r,, (x, i)) + F + ei(x, t), (3.2) 



where /'(x, r\\) = f\\ ( x, r\\ , rj_(x, r||)) and r^(x, r\\) is de- 
termined by Eq.(3.1). It is quite plausible that, after 
averaging over all f , the correlations in /' will also be 
short-ranged, albeit different from those of f, since the 
dissipative dynamics will avoid maxima of the random 
potential, effectively reducing the average forces. In that 
case, the equation reduces exactly to the model studied 
by NSTL and NF. Thus, the scaling of longitudinal fluc- 
tuations of the FL near threshold will not change upon 
taking into account transverse components, and the ex- 
ponent relations ( pr^ - |l.ll| ) hold for Model A as well. We 
expect this argument to hold even for Model B [Eqs.( p^ ] 
as long as V^J^i ^ ^x^|| ' '^^ when C±_ < C|| ■ 

For the interface model, it is possible to show that 
v{F) is a single valued function using the "no passing 
rule" of Middleton jj] . The rule states that no interface 
(or CDW) can overtake another, if initially every point on 
the flrst interface is behind the second one. This rule does 
not apply to the vector model: It is in principle possible 
to have coexistence of moving and stationery FLs, allow- 
ing for the possibility of a discontinuous (multi-valued) 
v(F). However, since a moving line samples an arbitrar- 
ily large region in the medium, it is plausible that the 
velocity self-averages at long times, resulting in a sin- 
gle valued v(F) (i.e., no hysteresis). However, flnite-size 
systems do suffer from such hysteresis which adversely 
affects numerical simulations of the model. These issues 



are further discussed in Sec. VII 



Several exponent identities can be deduced from the 
form of the linear response. 
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dra{q,uj) 



(3.3) 



in the (q, ~> (0 , 0) limit. Due to the statistical sym- 
metry of Eqs.( ^.lC| ) under the transformation r± — > ~r±, 
the linear response is diagonal. Let us first set a; = and 
examine the static response: An additional static force 
£(q) with zero spatial average (no q = component) can 
be exactly compensated by the coordinate change 

The distribution of f does not change in the primed co- 
ordinates. Thus, the static linear response has the form 



Xa7(q, W = 0) = da-, 



1 



Kq^ 



(3.4) 



Since e|| scales like the applied force, the form of the 
linear response at the correlation length ^ gives the ex- 
ponent identity 



C|| + i/i^ 



(3.5) 



Considering the transverse linear response seems to imply 
C_L — C|| ■ However, as will be shown below, the static part 
of the transverse linear response becomes irrelevant at the 
critical RG fixed point, since z± > Z||. This is consistent 
with the expectation that the dynamics is responsible 
for the distinction between longitudinal and transverse 
modes. 

Why are the relaxational dynamics different in the two 
fluctuation directions near depinning? The answer can be 
traced to a simple symmetry argument, which requires F 
and V to remain parallel, i.e. 



F(v) = F{v)i 



(3.6) 



where v = v/u, and F is some (scalar) function which 
depends on only the magnitude u, of velocity. For small 
deviations around v = ue||, we thus obtain (see Fig. ^ 



.^_dF 

9f II dv ' 
dF± _ F 
dv± V 



(3.7) 
(3.8) 



dF, 



V 



dv, 



dv^ 



dF 



1 



FIG. 2. A graphical demonstration of Eqs.(|3.7f- p.8| ). When 
a longitudinal force is applied, the direction is not changed 
and all changes are in the magnitude F{v). For a transverse 
force, F does not change to linear order in dF± , but v changes 
direction to remain parallel to F. 



These two derivatives clearly scale differently in the 
w — *■ limit, which causes a separation of relaxation time 
scales, as shown below. 

Now consider the response to a spatially uniform (q = 
0), but time-dependent, external force e{t). The leading 
term in the dynamic response is intricately connected to 
v(F): When a slowly varying uniform external force £{t) 
is applied, the FL responds as if the instantaneous ex- 
ternal force F -I- £ is a constant, i.e. it moves with the 
average velocity 



Ov 

(dtra) = Va{F + e) W Va{F) + -Q^^J- 

Therefore, near the depinning transition, 

X||(q ^0,uj)~ . I r^t 

—iu;(di'/av) + C>(uj'') 

1 

'iLu{Fc/v) + 0{uj^)' 



XJ.(q = 0,w) 



(3.9) 

(3.10) 
(3.11) 



Eq.( p.4[ ) can be combined with the above to yield a 
Taylor expansion of the inverse linear response around 
{q_,oj) ~ (0,0) that reads 

X|[\q, w) ~ Kq^ - iuj{dF/dv) + h.o.t., (3.12) 
X7^(q,w) ~ Kq^ - iuj{Fc/v) + h.o.t. (3.13) 

The zero of x^^ in the complex cj— plane for a given value 
of the wavevector q gives the relaxation time of the cor- 
responding mode. Thus, the relaxation times of fluctua- 
tions with wavelength ^ are 



, dv 
dF 



Fr. 



which in turn yield 



+ (2|| ~ 2)iy, 
z± = Z|| + 1/iy. 



(3.14) 
(3.15) 



(3.16) 
(3.17) 



Thus, z± > Z|| as noted earlier. This difference arises 
entirely from the different scaling properties of dv/dF 
[- {F ~ Fc)^-'^] and v/F [~ (F - F^ f] near the depin- 
ning transition, as noted earlier. 



B. Model B 

A similar linear response analysis can be made for the 
more general case of Model B. The leading contributions 
to the static and dynamic part of the inverse linear re- 
sponse are given by 



Xa7(q'^ = 0) = Ka-fq^, 

Xaji^l = 0,uj) = -iuj{dva/dF^y 



(3.18) 
(3.19) 
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The relation between the external force and the drift ve- 
locity can in general be written as 

F(v) F{v,0) [cos(t){v,0)v + sm(t){v,0)x x v] . (3.20) 

Both F and in general depend on the orientation of v, 
parametrized by an angle 9 in the yz-plane. Then, for 
small deviations around v = Tjeii, 



(3.21) 



where 

All = dy{Fcos(j)), 

Ml = (9i,(Fsin(^), 

All — F sin — 9e (F cos (/)) , 

A22 — F cos (f> + d0{F sin (p) . 

The scaling of diagonal elements in the linear response 
are t he same a s in Model A. Therefore, exponent identi- 
ties (3.16-3.17) hold in the more general case of Model B 
as well. 



suppressed henceforth. Since Z = 1 independent of the 
realization of randomness, response and correlation func- 
tions can also be generated using the disorder-averaged 
generating function Z = J [c?r] [c?f] (exp(iS')). For exam- 
ple, the two-point correlation function is given by 

(r„(x,t)r^(x',t')> = / [dr][df]r„(x,t)r^(x',t')(exp(5)). 



and the linear response is 



i j [dv] [df]ra(x, t)r^{yi' ,t'){e-x:p{S)) . 



In order proceed, we discretize in x-space: r(x, t) — > 
Vi[t). Introducing two conjugate fields Ri(t), Ri(i), Z 
can be rewritten as 



(4.4) 



Z = j [dR][dR]exp(S'), 

5 = 5]lnZ,(R„R,) 

3 

-I j di^R,(t)-J,r.i.R^.(t), 



(4.5) 



IV. MSR FORMALISM 

We use the formalism of MSR ||l^ to compute response 
and correlation func tions for the dynamical system de- 
scribed by Eqs. (2.9). After some rearranging, we obtain 



^drj^ = j A' J„,(x-x')r,(x',i) 

-r„(x, t) + (x, r(x, t)) +F^+ e„(x, t), (4.1) 
where the tensor J is given by its Fourier transform as 



Ka^q ■ Introducing an auxiliary field 



f(x, t), the generating functional Z is given by 



Z = J [dr][dr]J[r] exp(5). 



(4.2) 



where 



x' 7 



S^iJ d'^xdtfa{y^,t) [rjdtr-y - Ka-yVy 

-F„-/„(x,r(x,t))-£„(x,t)}. (4.3) 

Clearly, this coarse-grained continuum picture of the 
system breaks down at length scales shorter than the 
core radius of the FL. Therefore, there is a natural cutoff 



A in q-space for the functional integrals in Eq. (4.2). Z 
can be used to generate response and correlation func- 
tions of r, since integrating over f gives delta functions 
that impose the solution to the equation of motion (4.1). 
The Jacobian J\r\ fixes the renormalization of Z such 
that the delta functions integrate to unity, and will be 



where Zj ( Rj , Rj ) is given by 



Z,^/[.r,]M..](exp/..[.R,(.)-r.W (4.6) 

+lV,{t)-{7^dtV,{t)~B.,{t) 

+v,{t)-i,(v,{t))-¥-e,{t)] 



Note that this factorization of the disorder-dependent 
part of the action to local functionals Z j is possible only if 
the random forces are independent at each site j, as as- 



sumed in Eq. (1.3). Z j can be evaluated by an expansion 
around the saddle-point approximation. The integrand 
of the exponential is a maximum when, for each j, 

i 

-^J,r.i.RO + (r,)^p = 0, 

i 

which has a solution R" = 0, R" ~ vt for all j. Here, 
V is determined self-consistently as a function of F by 
requiring (rj)j^p = vi, where averages ((. . .)mf) are gen- 
erated from Zj evaluated at the saddle point, which is 
identical for each j: 

ZM. = /Mr,]Mf,](exp./dtf,W (4.7) 
• {vdtr,{t) - vt + r,{t) - f, {r,{t)) - F -£,(*)} V 
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Zmf can be identified as the MRS generating fun ction 
for a mean- field (MF) approximation to Eq. (4J), ob- 
tained by setting Jaj{^ ~ x') = 5a'yN~\ where = 
/ d'^x. (The first term in the RHS of (Q is then self- 
consistently equal to (r)MF = vi.) Redefining the field 
variables as R — > R + vt, iR —>■ R (for notational sim- 
plicity) , the expansion for In Zj I R^ , Rj ) is given by 



\nZj I Rj ,Rj 



ma'.ric 



a I Sq — 1 



^ / in n ^W^,a(iasJ n <<^J"(*-^ 

<V{™„.„„}({iasJ;{Q})- (4. 



The vertex functions V are obtained by evaluating deriva- 
tives of In Zj with respect to the fields at the saddle point, 
and are given precisely by connected correlation an d re - 
sponse functions of the MF system decribed by Eq. (4.7): 



V, 



{mQ,n„} {{^asa,}] {tas'^}^ 



(4.9) 



nn 



n n '^j"(^"''c.) 



MF,( 



Thus, once the mean-field system is solved, correlation 
functions of R, R can be studied through a momentum 
space RG treatment to obtain the scaling exponents of 
the fields in the long-time, large wavelength (hydrody- 
namic) limit. R and R are like coarse-grained forms of 
the original fields r and f since all correlation functions 
of r, f are equal to corresponding correlation functions 
of R, R in the hydrodynamic limit Therefore, it is 
sufficient to find the scaling behavior of R, R to deduce 
the desired critical exponents. 



V. MEAN FIELD THEORY 



In this section, we calculate response and correlation 
functions of the local system described by Zmf, which 
gives the vertex functions in the diagrammatic expansion 
of S. We will only need to calculate the leading terms as 
higher order vertices will turn out to be irrelevant in the 
critical region. Due to the averaging, Zmf is identical at 
all sites j, and it is sufficient to examine a single point. 
Setting r{t) = rj{t) — vt, and e{t) = £j{t), the equation 
of motion becomes 




(a) 




FIG. 3. (a) The effective potential Vaft- The random part 
(not shown) superimposed on the paraboloid slides with ve- 
locity —V. (b) A cross section of VcS- The particle stays in a 
local minimum P for a time of 0{v~^), after which the mini- 
mum disappears and the particle finds another local minimum 
P' within a finite time. Time averages are dominated by the 
slow portion of the motion as w ^ 0. 



V { ^ + Vc.) = -r^ + + (vi + r(t)) + e„(<). 



(5.1) 



F is determined as a function of v self-consistently by re- 
quiring that (f)MF — 0. The scaling behavior of Fmf(v) 
near threshold can be determined from the following ar- 
gument: For V <C r]~^F, the particle follows a local min- 
imum P of the effective potential 



14ff(r,t) = y(x,vi + r(t))-|- 



F.f(i). 



A representative snapshot of V^ff, which consists of a 
paraboloid centered at f = F with a superimposed ran- 
dom potential, is shown in Fig. |^. The position of the 
local minimum P shifts with a velocity of 0{v) as time 
progresses. Eventually, P disappears at a saddle point 
as it is pushed up the sides of the hyperparaboloid. At 
this moment, the particle quickly moves to a new local 
minimum P', after which it starts following the slow mo- 
tion of P', as shown in Fig. For scalloped random 
potentials with discontinuous derivatives at the saddle 
points, the particle starts moving with a velocity of 0(1) 
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(i.e., independent oi v a.s v ^ 0) as soon as P disap- 
pears, and reaches the vicinity of P' in 0(1) time, giving 
the result /3mf = 1- (In contrast, for smooth potentials, 
there is a w-dependent acceleration time just after P dis- 
appears, which contributes to the critical dynamics and 
gives /3mf = 3/2. We have also numerically inte- 

grated Eq. {bA) (for Model A) to verify that /3mf = 1- 
Next, we proceed to compute vertex functions 

"^'{mQ.n„} (^{^qsq}; {t'as' }) ™ ^^"^ pcrturbativc expansion 

of S*, which correspond to response and connected cor- 
relation functions of the MF theory, in increasing or- 
der in the field variables R, R. From now on, we set 

m = TO|| + mj_, and n — n\\ + nj_. 



A. Average position (m=l, n=0) 

By construction (f)[^p = 0, but we prefer to expand 
around the true F(v) instead of the mean- field value of 
the force Fmf(v). Since the effect of an additional uni- 
form static force F — Fmf(v) can be fully counteracted 
by a shift in r, this does not affect connected correlation 
or response functions. Thus, the only effect of this shift 
is to produce an additional term X^J^'' ~ Fmf(v)] • R^ 
in S, which only has a, q — component and does not 
directly enter the renormalization of higher order terms. 



B. Linear Response (m=l, n=l) 



5FMFa(v) 



+ 0(cj2). (5.4) 



Since /3mf — 1, the linear response tensor will have the 
form 



Xiuj) = l + iuj 



.^21 



(5.5) 



where Aa-y approach constants as u ^ (cf. Eq.(3.21)). 
For Model A, x{uj) is diagonal due to symmetry, and 
Ai2 = A21 = 0. 



C. Nonlinear response (m=l, n>l) 

Assuming that Fmf (v) has a Taylor expansion around 
V = we|| for u > 0, we can expand the RHS of Eq. ( |5.3| ) 
to obtain the nonlinear response of the model. The lead- 
ing term in the low- frequency limit is proportional to a;", 
and it is straightforward to show that the contribution of 
these terms to S is 



1 



nil \n±l 



(Txdt 



9"F„(v) 



9"iii;||9"it;_L_ 
xi?„(x,t)(ati?||)"ii(9ti?i)^ 



(5.6) 



These terms are irrelevant at the RG fixed point, as we 
shall show later. 



The linear response is given by the rank 2 tensor. 

We are only interested in the low-frequency form of the 
Fourier transformed linear response x(w), i.e. when e 
is slowly varying in time. In this case, we can write 
£(t) ~ £0 + neglecting terms proportional to £. To 
find the response ?£ [t) , let us define 

Y'g-{t) = ve{t) - £0 - eit - Fmf(v) + Fmf(v + £1). 
Taking a time derivative and using Eq. (|5.lD , we obtain 

/?(?£ + V + £1) = -fg + Fmf(v + £1) 

-f f((v + £i)t + r^ -Fe), (5.2) 

where Fg = Fmf(v + £1) — Fmf(v) — £o- But now, 
(r^) = by definition of Fmf- (The random force f is 
evaluated at points shifted by a constant amount Fe, but 
this has no significance upon averaging over randomness.) 
This gives 

{re{t)) - e{t) + Fmf(v) - Fmf (v + e{t)) + 0{i). (5.3) 
Expanding Fmf(v + £) for small £, we obtain 



D. Two-point Correlation Functions (m=2, n>0) 

At low velocities, the particle spends most of the time 
near a local minimum, jumping abruptly to the next 
one when this minimum disappears. Therefore, the time 
scale associated with the correlation functions is given by 
the temporal separation between two consecutive jumps, 
which scales as 1/v. In the w — » limit, the correlation 
functions depend on t only through the rescaled time 
variable u = vt, since the positions of successive min- 
ima near threshold are determined by energetic consid- 
erations, and do not depend on v. (The correlation func- 
tions may also depend on the drift direction v. We shall 
suppress this dependence for notational brevity.) Let us 
define 

{r^{t)r,{t'))^^^^^Cc,,{v{t-t')). (5.7) 

Since successive positions of the local minima are uncor- 
related, we expect that Ca-y(u) decay quickly as a func- 
tion of u = vt for \u\ > 1. By definition, 

^C^{u) = Ci2{u)^C2i{-u). 
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As a result of the abrupt jumps from one minimum to an- 
other, Ca^{vt) have a discontinuous derivative at the ori- 
gin, rounded at a scale of 0{v). In Model A, Cx (w) = 
due to symmetry. 

The only other important terms in the effective action 
S involve the series m = 2,n = n\\ > 0,n± — 0. All 
vertex functions associated with this series are given by 
the response of connected correlation functions to longi- 
tudinal forces. These response functions are intimately 
related to the two-point correlation hmctions Cq.^(u) by 
the following argument: Static forces only change linear 
response, and do not affect connected correlation func- 
tions. For a slowly varying external force e(i)e|| , however, 
the system will respond as if the instantancus velocity is 
{v + e)e||. Neglecting terms proportional to e. 



MF,, 



C^^{{v + e){t-t')) + 0{e) 
Ca-y {v(t - t') + e{t) - e{t')) . 



Now, Taylor expanding Ca-y around v{t — t') and taking 
successive functional derivatives with respect to e, we fi- 
nally obtain the contribution of this series to S as 

^ = f d'^xdtdt' Ra{x,t)Ry{x,t') (5.8) 

n=l 

X {/„^,„ {v{t - t')) (X, t) - (X, t')] " , 

where Uaj.niu) is the nth derivative of Cajiu). 

The vertices with m = 2, n± > and m > 3 are all 
irrelevant, as shown in the next section. 



VI. SCALING AND RG 



The terms in S that are up to second order in the fields are 



So = - J dt d'^x [F - Fmf(v)] • R(x, t) 

-CM j-i(q)-xH 



1 
2 



q,cj 



R(-q,-c^)' 
R(-q,-^) 

S, 







R(q,w)' 
R(q,c^) 



(6.1) 



where J^^{q) = (^Saj — Kajq^) ^ ~ Saj + Kajq^ for small q. For notational brevity, we use ^ to denote / ^^■J)d 
Using Eq.(5.5), the quadratic form in the action can be written as 



1 

2 



i?l|(-q, ~uj) 
i?i(-q, -uj) 

^11 (-q, 



■ Q(q,^) 



^l|(q,u;) 

^±(q,^) 

i?l|(q,w) 



where 



Q(q,'^) = 



-~C\\{uj) -iCx(tj) Kiiq^-iuAii vKi2q'^ + iuAu 

-^Cyi-ij) -C^{lo) 
Kiiq^ + iujAii K2iq^ + iw^2i 



K2iq^ - iujA2i vK22q^ - iL0A22 



vKi2q^ — iujAi2 vK22q^ + iw^22 













(6.2) 



Neglecting all higher order terms in the action, we 
arrive at a Gaussian theory, in which different Fourier 
modes are decoupled, an d w hich can be solved by in- 
verting the matri x in Eq.( |6.2| ). (See Appendix ^) The 
quadratic action (BT) remains invariant under the scale 
transformation 



fox, 



t b'^t, 



R\\ -> 6-2-d/2^^| ^ ^ j-2-d/2^^^ 

y^b-'^/\, F-FuF^b-^/^iF-FuF): 

(6.3) 

except for terms proportional to Ari2 and ^22 which 
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vanish at the depinning transition as v — > 0+. For 
d > 4, ah higher order terms in S decay away upon 
rescahng, and we recover an asymptotically quadratic 
theory with critical exponents P = 1, = 2, v = 
2/d, C|| = (4 - d)/2, C± = 2 - d. The remaining ex- 
ponent, z±, can be found by comparing the static and 
dynamic parts of the transverse linear response. This 
gives z± — 2 + d/2 = zn + l/v, as shown previously by 
the exponent identity (1.13). The exponents related to 
longitudinal fluctuations, not surprisingly, are identical 
to corresponding exponents in the interface problem 
However, we have also calculated new exponents char- 
acterizing transverse fluctuations. We see that even the 
simple Gaussian theory exhibits anisotropic exponents. 
At d = dc = 4 dimensions, the scaling dimension of 
changes sign and we cannot neglect its higher powers 
anymore. Simple dimensional analysis indicates that the 
only higher order terms in S which become margin al a t 
d 



dc involve vertex functions Uaj.n, given in Eq.(5 



This series can be summed up over n, together with the 
n = term included in the Gaussian theory, to yield 

d'^xdtdt'Ra{x,t)R^{K, t') 

X {v{t-t')+ R\\ (x, t)-R\\ (x, t')). (6.4) 

All higher order terms in S are formally irrelevant since 
they involve additional powers of i?^, or R^, whose 
scaling exponents are less than zero. 

For d < dc, the vertex functions Ua-f^n become more 
and more relevant for increasing n under the rescaling 
( |6.3[ ), and the fixed point moves away from the Gaussian 
theory. In d = 4 — e dimensions, we look for new fixed 
points with different scaling properties: 



R\\ ^ b'^w R\\, RjL^bi^R^, 
i?ll ^ 6^i|-''i?l|, RjL^b'^^-'^R±, 
F- Fmf ^b-^/^iF- Fmf), w^b-f^l^'w. 



(6.5) 



To calculate the new exponents to first order in e, we 
employ a one-loop momentum shell RG scheme, tr eatin g 
all non-Gaussian terms in the action (i.e. U in Eq.(5.8)), 
as a perturbation. Perturbative calculations proceed by 
expanding (e'^)o, where (• • ■)q denotes averaging with re- 
spect to the Gaussian action 5*0, in powers of U. A renor- 
malization transformation is then constructed as follows: 
(1) Perform the averages only over short wavelength fluc- 
tuations R^,R^ with wavenumbers A/6 < |q| < A, 
where b — e*^. The resulting coarse grained action is 
perturbatively given by 



5'^ — 



mi 



\m 



> 

0,c 



0(U^). 



(6.6) 



(2) Apply the rescaling transformations given in ( |6.5| ), 
bringing back the short-distance cutoff A to its original 
value. (3) The exponents are then determined from the 
fixed points associated with the RG flows of the the ac- 
tion. Since Models A and B are characterized by distinct 
fixed points, we shall discuss them separately. 



A. Model A 



In the low- frequency, small- wavevector limit, the effective action for Model A is 
5^^) = - J dtd'^x[F ~ FMFiv)]R\ii^,t) 

{i?||(-q,-c^)i?||(q,c^)(ifg2 _^^An) + i?±(-q, -Cc;)i?^(q, c^)(i^g2 _ .^^^^/w)} 



+ i y" d'^x dt dt 'i?|| (x, i)i?|| (x, t')C\i {v{t - t') + (x, t) - (x, t')) 
+ i f d'^xdtdt'R±{x,t)R±{x,t')C±{v{t-t') + Rn{x,t)~ R\i{x,t')) . 



(6.7) 



The Gaussian part has the correlation functions. 



{R\\{-(l,-i^)R\\{q,uj))a = — — , 

1 

(i?±(-q,-cj)i?i(q,tj))o = 



AV - iuiA22/v' 



(6.8a) 
(6.8b) 



(-Rii(-q,-'^)-R|i(q,^^))o = 
(i?_L(-q,-t^)i?±(q,w)>o = 



xV + (^^ii)2' 

K^q^ + iLuA22/vy 



(6.8c) 
(6.8d) 



The vertex functions Uaj^n = for a 7^ 7, and these 
terms are not generated by the RG transformation. The 
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renormalization of remaining vertex functions C/||^„, and 
U±^n for n > can be recast into a functional renormal- 
ization of C\\{vt) and C±{vt), provided that vt and R\\ 
scale in the same way, i.e. C|| — zy — P/iy. This rela- 



tion can be independently obtained from Eqs. (3.5) and 
( |3.16 ), derived in Sec. II] from more general (and nonper- 
turbative) arguments. The renormalized vertex functions 
are then obtained from successive derivatives of C{vt) as 



C/a,«(") = C(")(u). 



(6.9) 



This ensures that the form of Eq.(6.7) is retained un- 
der renor maliza tion, albeit with renormalized parame- 
ters. Eqs. ( 6.8c ) and ( 6.8d ) suggest that Ca{vt) may be 
interpreted as temporal correlation functions of an effec- 
tive force generated by the quenched disorder. 



The renormalization of some terms in Eq.(3/7) do not 
get any contribution from the momentum shell averaging 
step, giving rise to additional exponent relations that are 
correct to all orders in the e expansion. The first relation 
is due to the fact that F never appears explicitly in any of 
the contractions or higher order vertex functions. Thus, 
the renormalization of the term proportional to _F — Fmf 
can be written as 



djF - Fmf) 
d£ 



iz«+en){F-F, 



MF 



constant, (6.10) 



where "constant" refers to an expression that does not 
involve F. This RG flow equation can be rewritten as 



d{F - fc) 
d£ 



(z||+f?[|)(^^-F,), 



(6.11) 



with a suitable choice of Fc. Hence, higher order correc- 
tions may shift the threshold force, but do not influence 
the scaling of F — fc- This implies that 



9|| - l/iy = 0. 



(6.12) 



Furthermore, there are no contractions that contribute 
to the renormalization of K or A22. Thus, 



'II 



+ 2|| +c 



0, 



9^+C±+ (ijv = 0, 



(6.13) 
(6.14) 



respectively. As a result, all critical exponents are deter- 
mined in terms of C|| , C± ^-nd zy . These exponents can 
be computed by constructing RG flow equations for the 
remaining parameters. 



1. Renormalization of Ca 



dC»{u) 



=[e + 2^^ll+2(z||-2)]q|(ii)+Cpq(7/) (6.15) 

-if,{[qi(u)]^ + [q|(«)-q|(o)]c[|'(^)}. 



dC^{u) 
d£ 



[e + 26l_L + 2(z|| -2)]C_l(u) 



(6.16) 



+ Cii u {u) - Ka { [q, {u) - q, mc'i {u) } . 

The constant Kd = SdA'^^'^ /[i27rYK^], where Sd is the 
total solid angle in d-dimensions. Primes denote deriva- 
tives with respect to u. Terms proportional to uC'^{u) 
arise from the rescaling of u. We look for fixed-point so- 
lutions C*(u) that decay to when |m| is large, since they 
are related to correlation functions of the system, which 
are expected to vanish for large time differences. 

Not surprisingly, the functional recursion relation for 
C||(u) is identical to the one obtained in Rcf. In fact, 
all higher loop corrections are identical as well. This 
is in excellent harmony with the argument presented in 
Sec. Ill, and allows us to use the resu lts of NF. Setting 
dC^^/dl — 0, and integrating Eq. ( 3.15| ) from u = —00 to 
00, we get 



[6 4-2011 +2(z|| 



2)-C||] j q|(u)dw = 0. (6.17) 



Provided that the RG flows go to a fixed-point solution 
with / C* 7^ 0, this implies that C|| = t-2[2- {z\\+e\\)]. 
The mean-field correlation function satisfies this integral 
condition for both random-field and random-bond dis- 
order, since C is essentially insentitive to the value of 
the random potential between consecutive local minima 
P, where the line moves quickly. There are other fixed 
points with / C* =0, but they are irrelevant for our 
discussion. Thus, from Eqs. (3.5) and ( 6.12| ), we obtain 



C|| = e/3, 
3 



6 - e 



(6.18) 
(6.19) 



NF prove that these results are correct to all orders in e, 
by showing that the contributions to the renormalization 
of C|| (u) from higher-order terms is a complete derivative 
with respect to u. Upon int egrat ion over u, such higher 
order terms do not alter Eq.( j.fij ), leaving the exponents 
unchanged. 

Using C|| — e/3, an implicit solution for C*^{u) is ob- 
tained as 



qT(u)-qT(o)-qT(o)in 



cm ] ' 2 



After performing the momentum shell integration and 
rescaling, details of which are given in Appendix we ar- 
rive at the recursion relations for the renormalized func- 
tions Ca{u): 



where uq = y/3KdC*{0)/e. q*(0) is arbitrary, and can 
be changed by a rescaling of the fields . Expanding 
the logarithm for small m, we see that there is a kink at 
the origin, as 
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q(o) 



lul 1 / u 

1- — + :t - 

uo 3 \uo 



+ 0(|u/uo|3). (6.20) 



For \u\ 3> uq, the fixed point solution behaves hke a 
Gaussian, and 



C||(w) « C||(0)exp 



2^0 



We next examine the fixed-point solution C^(u), 
which is the new element of our computation. Setting 
dC^/di — and looking at the limit u ^ 0+, we get 
(assuming that Cl{0+) ^ 0) 



[£ + 261^ + 2(z|| -2)]C1(0+) = 



Combined with Eqs.(U), ( ^.16D , and ( |6.14| ), this result 
implies 



5e 



(6.22) 



In Appendix |c[ we show that this result is in fact cor- 
rect to all orders in e since there are no contributions 
to C^(0+) from momentum-shell integration. The fixed 
point solution (for m > 0) satisfies the equation 



^ln|CI(.)| = 4 



du 



cm 



(6.23) 



3.21) Upon integrating twice, Eq. ( 6.23| ) leads to 



CUu) 



-Cl'(0+) J du' exp i --^ y" du" — 



q(^.")/q(0) 



(6.24) 



where (0+) is arbitrary in the same sense as C||(0). 
For |w| > Mo, Eq.(|6.23|) gives 



Cliu) 



C-^exp 

u 



2ul 



(6.25) 



where C is a constant related to (0+). The numerical 
solution for the fixed point functions C* (u) are shown in 
Fig. ^. The qualitative features of C|* and are sim- 
ilar: both have a discontinuous derivative at the origin, 
and decay as a Gaussian for large values of However, 
note that their scaling dimensions differ by . 




FIG. 4. Fixed point functions C|*(m) (solid line) and C'^{u) 
(dotted line) , normalized to yield 1 at the origin. Their values 
for M < (not shown) are found from C^{u) — C^i—u). 

The exponent ^|| = e/3 can also be obtained by naive 
dimensional arguments: In dimensions d > 4, the ran- 
dom force can be expanded as /|| (x, r|| ,r±) = /|| (x, 0, 0)-|- 
0(r||,rj^). Since both r|| and rj_ have negative scal- 
ing dimensions (C||,C-L < 0), the correction terms can 
be ignored. The random force scales as b"'^^^ under a 



scaling x — )■ 6x, leading to the Gaussian roughness of 
C|| = 2 — d/2. A similar scaling argument applied to 
Eq.(p]) leads to C-L = C|| - d/2 = 2 - d. For d < 4, the 
scaling dimension of r|| is positive, and higher powers of 
r|| in an expansion of /y (x, ry , r±) are more relevant. It is 
then reasonable to assume that in this case the statistical 
properties of /|| at large r|| arc crucial. If uncorrelated at 
large separation, /||(x, r|[,0) scales as fo^^'^+^n'/^. When 
equated to 6''ii~^ for the scahng of V^ry, this leads to 
^y = e/3 in agreement with the RG treatment. Essen- 
tially, the statement regarding the non-renormalization 
of J duC|[(u) justifies the above "naive" scaling. How- 
ever, a similar reasoning from Eq.( [3.l| ) would have con- 
cluded C-L — C|| = —{d+ C||)/2, in disagreement with 



Eq.(3.22). In this case, J duC±{u) is renormalized, but 
C±{0) is not; suggesting that despite the presence of rele- 
vant higher order powers in the expansion of /^(x, r|| , r^) 
around r = 0, the scaling properties are still controlled 
by /i(x, 0,0). We have no physical motivation for this 
rather curious conclusion. 



2. Propagator Renormalization 

The only one remaining exponent is z\\, which can be 
obtained by examining the renormalization of An. One- 
loop contributions arise from the n = 2 term in {U)q, 
which is 

oo oo 



iy'd'^x J dt J dt'7?y(x,i)7?y(x,t') 



-oo — oo 



X[i?y(x,t)~i?y(x,t')]2q'(,;(t-t'))- 

Replacing [i?y(x, t')]^ with [i?y(x, t)]^ does not change 
the integral. Thus, upon further manipulation, this term 
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in the action can be written as 



d'^x dt dt'7?||(x,t)i?||(x,i')^||(x,t) 



-OO — OO 



X (x, t) - i?|| (x, <')]C||' {v{t - t')) . 

Since a contraction forces t and t' to be within 0(1) of 
each other, and we are only interested in the first time 
derivative, we can substitute i?||(x, t) — _R||(x, t') « (t — 

t')dtR\\{K,t). Now, contracting i?||(x, i') with i?||(x, t) 
and integrating over the momentum sheU, we obtain a 
contribution to An equal to 



se 



(27r)''Ai 



dt t e' 



-Kts?tlAx\(~iii 



"qi'(ft). (6.26) 



The minus sign comes from the opposite overall signs of 
m— \ and m — 2 terms in Eq. (3.7). For w — > 0, we can 
set the argument of Cy' to zero. However, this causes a 
problem: Cy' has a term proportional to 5{vt) in the low- 
frequency analysis, this term diverges as \/v for — > 0. 
This apparent divergence cannot be avoided within the 
low- frequency analysis we have used so far. The propa- 
gator is sensitive to high-frequency behavior of the vertex 
functions. Careful analysis of the high frequency struc- 
ture of C|' shows that the terms that contribute to the 
diverging part of C|'(0) do not enter the renormalization 
of the propagator. (See Appendix ^.) This is essentially 
due to the causal nature of the response: Perturbations 
right after a jump do not influence the motion before the 
jump. The correct way to avoid these divergent terms 
within the low- frequency analysis is to use C|''(0^) in- 
stead of C"(0). Near the fixed point, this can be calcu- 
lated to 0(e) from Eq. ( |6^ ) as O||"(0+) = 2e/(9Xd), 
resulting in 

v4< = All - 5iAiiKdC{ {0+) = Aii[l - <5f(2e/9)]. 

Finally, after performing the rescaling, we obtain the re- 
cursion relation 



dAi 



de 



- ^ An[0» + G - 2e/% 



which yields 



zn =2-2e/9 + 0(e2 



(6.27) 



(6.28) 



B. Model B 



The presence of off-diagonal terms in the action 
changes the critical scaling properties of Model B. The 
nonzero contractions that appear in the momentum shell 
integration in this case are (cf. Appendix ^) 



(i?l|(-q,-L^)i?||(q,c^))o = ^^^^ : , (6.29a) 

(i?i(-q,-c^)i?||(q,c^))o-— , (6.29b) 

(i?,l(-q,-c.)i?||(q,c.))o = -^^2^^^^, (6.29c) 



where 



K = A12/A22, 

K\\ = Kii 
p\\ = All + KA21, 
C{uj) = C|| {uj) + KRe[Ox (uj)] + K^C±{uj). 



In addition to the nonrenormalization relations ( 6.12 



6.14), the nonrenormalization of K21 or A12 dictates that 
6*11 = e±. (6.30) 
This immediately implies the exponent identity 

a = C[i - 1/^- (6.31) 



The naive scaling argument based on Eq.( |3.l| ) gives 
an equivalent result when the scaling dimension of 
dr^/dr\\ (Cj^ — ^||) is equated to the scaling dimension 
of /_i(x, r|| , 0) [—{d + Cy)/2]. The naive argument works 
this time, since J duC±{u) remains finite at the fixed 
point (see below). 

Under this rescaling, k and i^n remain unrcnormalizcd, 

and the renormalizations of p|| and C determine the re- 
maining exponents C|| and zy . The recursion relations of 
vertex functions Ca are more complicated, but there is a 
relatively simple fixed point solution with 

C*{u) = 4C'|]'(u) = 2kO* (it) = Ak^CHu). (6.32) 

Furthermore, C{u) satisfies a rec ursion relation identical 
to that of 0|[ (u) given in Eq.(|T|). This result once more 
shows that longitudinal fluctuations, whose correlations 
are given by Eq.( |6.29c ), are not altered by the introduc- 
tion of transverse fluctuations even in the more general 
case of Model B. 

The renormalization of p|| also gives results very simi- 
lar to that of Model A, with the substitutions O" C" 
and ^11 — !■ Thus, the RG analysis gives the same 
exponents C|| — e/3 and Z[| — 2 ~ 2e/9 -f O(e^). Further 
details appear in Appendix 

If the Hall Angle is sufficiently small, the FL can not 
distinguish between zero and nonzero angles. Therefore, 
the effective roughness and dynamic exponents at small 
length and time scales should be given by the Model A 
fixed point. Note that n = tan in a n isotropic system 
with nonzero Hall angle (cf. Eq.(3.21)), and k is in gen- 
eral strongly related to the macroscopic Hall angle. Thus, 
K ^ 1 when the system is almost Model A-like, and its 
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nonrenormalization determines the cross-over behavior to 
the Model B fixed point: Under renormalization with 
Model A exponents, the system remains near the Model 
A fixed point until the ratio C±/C\\ increases to 0(k~^), 
as the Model B fixed point is approached. Isotropic ef- 
fective exponents appear in this crossover regime. The 
length scale at which the behavior crosses over to the 
Model B is roughly given by 



(with Model A exponents for C,a), i.e. the anisotropy is 
noticeable when the angular spread in the direction of a 
typical avalanche is of the order of (j). Thus, for the FL, 



6 



The range tq, of the impurity potential was kept con- 
stant. This construction creates a random scalloped po- 
tential landscape, eliminating any additional crossover 
effects that could arise from a smooth potential. 

Unless noted otherwise, all presented results were ob- 
tained using a grid size Aa; = 1, and a time step 
Ai = 0.02, in order to optimize computational con- 
straints. (Smaller values of Aa: or Ai did not lead to sig- 
nificant improvements.) Free boundary conditions were 
preferred over periodic ones since scaling was observed 
over a larger range of length scales in the former case. 
Other simulation parameters were K ^ \, w = \, ~ 1, 
Smax = 2. We expected a threshold force close to 1 for 
these parameters. A summary of our findings is presented 
below. 



which diverges as </> ^ 0. When ^ < , the anisotropic 
fixed point is never approached. Thus, the true critical 
region can be very small and difficult to observe for small 
Hall angle. 



VII. NUMERICAL WORK 

In this section, we present and discu ss th e results ob- 
tained by numerically integrating Eqs. (1^), providing a 
test of the analytical results presented so far. There are 
several difficulties associated with numerically studying 
critical behavior in a finite system slightly above thresh- 
old. In order to obtain meaningful statistical averages 
one must wait for the system to reach a stationary state. 
However, for any reasonably broad distribution of pin- 
ning forces, the system always gets pinned after a time 
~ e 



{F-Fa) L ^ where L is the linear extension. Therefore, 
in order to probe the critical region, it is necessary to go 
to very large system sizes. 

The necessity of integrating big systems, and the large 
computational cost of implementing quenched disorder, 
forced us to restrict numerical simulations to d = 1, in 
any case the most physically relevant dimensionality. We 
were further motivated by the expectation that some ex- 
ponents were calculated to all orders in e, and thus could 
be checked even at e = 3. 

Integrations were carried out as follows: Coordinates 
X and t were discretized, but the position r was left con- 
tinuous. For each x, the value of the random potential at 
point r was determined from a superposition of arttactive 
impurity potentials 

U,[r')^\s,{r'^~rl)Q{r,-r% 

where Q is the step function and r' is the distance from 
the center of the impurity. The impurity centers were 
randomly placed with a density w] their strengths Si were 
randomly drawn from a uniform distribution [0, Smax)- 
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FIG. 5. A plot of average velocity versus external force for 
a system of size 2048. Statistical errors are smaller than sym- 
bol sizes. Both fits have three adjustable parameters: The 
threshold force, the exponent, and an overall multiplicative 
constant. 

The velocity exponent /3 can be extracted from a plot 
of velocity versus external force. Such a plot is given in 
Fig. H for a system of size L = 2048. Each data point was 
obtained by a time average over 10^ time units and took 
about 30 hours of CPU time on a Silicon Graphics R4000 
workstation. The best power law fit gives an exponent 
/3 K. 0.31, but a weaker logarithmic dependence, which 
corresponds to /? = 0, seems to provide a better fit to the 
data. The conclusion is that higher order terms in v give 
very large corrections to the scaling of u, since either /3 
is very small or exactly zero. /? = would imply that 
Z|| — 1, a possibility discussed by NF for interfaces in 
1 4- 1 dimensions 
0.93 and 0.94. 



The threshold force Fc, is between 



The roughness exponents C|| , C-i 
equal-time correlation functions 



([r,(x,t)-r,(x',t)]2) 



are extracted from 
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10° 10' 10' 10^ 

Ix-x'l 

FIG. 6. A plot of equal time correlation functions versus 
separation, for a system of size 2048 at F = 0.95. The ob- 
served roughness exponents are close to the theoretical pre- 
dictions of (^11 =1, (^x = 0.5, which are shown as solid lines 
for comparison. 

Results for a system of size 2048, at a driving force of 
0.95 [{F-F^)/F^ « 10-2], are shown in Fig. |. The aver- 
ages were taken over a time interval of 10^, after waiting 
for all correlations to reach steady-state. The results are 
in overall agreement with the predicted values of the ex- 
ponents, even at e = 3. The slightly smaller value of C|| 
is expected, since determination of the roughness expo- 
nent from equal-time correlations becomes unreliable as 
the exponent approaches unity, and is inappropriate if it 
exceeds 1 . The deviations of transverse correlations 
from the scaling form are likely to be due to crossover ef- 
fects: The analysis of transverse fluctuations in the crit- 
ical region is correct only when v/F <C 1, because then 
the static part of the transverse propagator can be ne- 
glected. However, in our simulations v/F 0.4, suggest- 
ing that the critical region is much smaller for transverse 
fluctuations compared to longitudinal ones. 

In order to obtain an independent estimate of the dy- 
namical exponent Z|| , we also examined fluctuations in 
the spatially averaged velocity as a function of time. The 
resulting measurements were related to the previously de- 
fined exponents by the following argument Slightly 
above threshold, the motion of the line can be thought as 
a superposition of avalanches of various sizes Z, with an 
average lifetime Z^n and moment . Such avalanches 
occur if a portion of the line finds a region of size Z'^+^n 
with weak impurities. Thus, ignoring all power-law pref- 
actors, the probability of such an avalanche for / ^ ^ 
is 

P(0^exp{~G/e)'+^"}. 

Velocities at two separate times are correlated if there is 
an avalanche that is active at both times. Therefore, it is 
reasonable to assume that at large times, the contribution 
of an avalanche of size / to {v{t)v{0))^ is proportional to 



e"*/' " , once again neglecting power-law prefactors that 
depend, for example, on the typical number of active sites 
at a given time during the avalanche. The total contribu- 
tion of all avalanches is given by an integral over all sizes 
I with the probability measure P{1) dl. The leading-order 
time dependence of the exponent can be determined by 
a saddle point evaluation of the integral, resulting in 
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FIG. 7. Velocity correlations versus time, for the same sys- 
tem in Fig. 1^. A stretched exponential is a good fit to the 
data. 

where j = {d + ^||)/(z|| + d + (\\), suggesting a stretched 
exponential. The numerical results and the fit to a 
stretched exponential are shown in Fig. |^. (It should 
be noted that a comparable fit can also be achieved by 
a sum of two exponentials.) Assuming that C|| = Ij we 
arrive at z w 1.3, which is consistent with the value of 
(3 « 0.31 found from the velocity-force relation. Unfortu- 
nately, the data becomes noisy at larger values of t, due 
to the finite size of the time window used to extract the 
correlation function. The small value of r(« 10) makes 
it hard to predict the reliability of this estimate, since 
the power-law prefactors may be large and nonnegligible 
for such moderate values of t. Unfortunately, improving 
on this simple estimate is difficult as the determination 
of power-law prefactors requires a number of additional 
assumptions that are hard to test. Nevertheless, based 
on the accumulated numerical evidence it can be reason- 
ably argued that ^|| is between 1 and 4/3, the 0(e) RG 
prediction. 

Computed longitudinal exponents are also in good 
agreement with results from 1-1-1 dimensional interface 
depinning models. Numerical integration of Eq.( ^.lC| a) 
for an elastic interface (no transverse component) 
has yielded critical exponents ^ = 0.97 ± 0.05 and 
ly = 1.05 ± 0.1. Similarly, the force vs. velocity data 
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has been adequately described by both a velocity ex- 
ponent /3 = 0.24 ± 0.1 and a logarithmic dependence 
V ^ l/ln(F — Fc), which corresponds to /3 = 0. These 
results provide strong support for our prediction that lon- 
gitudinal exponents are unchanged when transverse fluc- 
tuations are introduced. However, it should also be noted 
that experiments and various discrete models of interface 
growth have resulted in scaling behaviors that differ from 
system to system. A number of different experiments on 
fluid invasion in porous media ||30|| give roughness expo- 
nents of around 0.8, while imbibition experiments [ pT|j32) 
have resulted in C « 0.6. Some of these results can be 
explained by the effect of anisotropy, which will be dis- 
cussed in the next Section. On the other hand, a discrete 
model studied by Leschhorn gives a roughness expo- 
nent of 1 .25 at threshold. Since the expansion leading to 
Eqs.(L2) breaks down when ( approaches one, it is not 
clear how to reconcile the results of Leschhhorn's numer- 
ical work with the coarse-grained description of the 
RG calculation, especially since any model with C > 1 
cannot have a coarse grained description based on gradi- 
ent expansions. 



VIII. DISCUSSION AND CONCLUSIONS 

In order to put the results we have found so far in 
better perspective, it is useful to discuss the effect of 
nonlinear terms that were ignored earlier, aspects of uni- 
versality, and possible generalizations to other systems. 
These issues are discussed below. 



A. Nonlinear Terms 



The leading order nonlinearities in Eq.(2.6) can be ex- 
amined by a gradient expansion, being careful to treat 
terms of 0{{dxr)^, {dxr)^dtr) accurately. After some re- 
arrangement, we arrive at 

y/1 + s||2 " ^ 2 " 2 

+AixS||S± +i^ + /||(x,r,S||,s^), (8.1a) 

r]dtr± 



+A2xS||S^ +/_i(a;,r, S||,s_l), (8.1b) 
where s\\ = dxr\\, s± = dxr±, and the random forces are 

(/||-S||/a:) , 



/II 



- [1 + V2 - SiV2]/± + s^L 
H 7^ sm ( 



/ (/-I - S±fx) 

f 1 = , COS ( 



[l-,s„V2 + ,s^V2]/||-.||,sx/^-,S||/. 
The remaining parameters are given by 



sm ( 



F 




Ai|| 


= -Fsin^c/), 


Ai± 


= -Fcos^ (j), 


Aix 


= -i^sin^^, 


A211 


= F sin (j) cos (j) 


A2± 


= — _Fsin0cos 


A2X 


= F cos^ 0. 



These equations of motion, and their generalizations to 

X G have thus been complicated by two factors: 

There are orientation-dependent terms, and the mean 

- 2 

square of the random forces = {fa ) also depend 
on the local orientation of the FL. By naive dimensional 
counting, it can be immediately seen that Ai|| and A211 
are relevant with respect to the fixed points we have dis- 
cussed for d < 4. In t he c ase of Model A (isotropic 
disorder with 



ridtrn 



0), Eq.(BJ^) further simplifies to 

f\\ - fxS\\ 



F 



F 



2 - 
Fs\\si_ + 



;.2a) 
.2b) 



Note that the two relevant nonlinearities vanish, and that 
Aa does not depend on orientation up to and including 
O(s^). Dimensional counting suggests that the remain- 
ing nonlinear terms are irrelevant and Model A exponents 
are valid for d > 1. Many more nonlinear terms become 
marginal d — 1, and the gradient expansion breaks 
down. It is unlikely for the critical exponents to change 
their value discontinuously at d = 1, although logarith- 
mic corrections to scaling exponents arc quite possible. 

The fixed point investigated here is unstable and only 
approached at the dcpinning force. Away from the 
threshold, critical scaling laws are observed at scales 
smaller than the correlation length scale ^. Beyond this 
critical regime, the behavior of Eq. ( ^.1] ) is similar to 
regular diffusion with white noise (a multicomponent 
Edwards- Wilkinson (EW) equation |Q), or the gener- 
alized KPZ equation (l|Jl|,^). A nonzero Ai|| of 0{v) 
is generated kinetically in this regime even if the system 
is initially isotropic with = 0, due to the terms on the 
left-hand side of Eq.(^.l|a). For d < 2, this nonlinearity 
is relevant, while for d > 2, a critical value Ac separates 
a weak-coupling region described by the EW equation 
from a strong-coupli ng reg ion described by the (general- 
ized) KPZ equation |l^ , [15|j3q ] . 

When 7^ 0, even in a fully isotropic medium, the rel- 
evant nonlinearities are nonzero, and the system is driven 
away from the "linear" fixed points. We discuss this and 
other possibilities next. 
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B. Anisotropy and Universality 

We noted earlier that anisotropy plays an important 
role in determining scaling properties near depinning, 
even in the absence of nonlinear terms. To fully under- 
stand the effects of anisotropy, including nonlinear terms, 
let us start by considering the simplest prototype of a 
FL oriented along the c— axis of a High superconduct- 
ing single crystal, such as YBCO. For simplicity, assume 
that the system is completely isotropic in the y—z plane, 
with = 0. Then, the motion of the FL is governed by 
Eqs.(^.2[), and the only important source of anisotropy 
is due to (/||^) = (fx)- This causes the mean 

square magnitude of /|| to depend on the local orientation 
as. 



A 



(/||'> + ((/.^>-(/||'» 



'II 



For interfaces, the depinning force is known to scale with 
the strength of the disorder i.e. - A^/f^-d). 

Thus, A II creates an orientation dependent depinning 
force il. 



2/(4-d) 
II 



Fr 1 



2 {m~{f\i') 



(8.3) 



This leads to a no nzero Aiy when the nonlinear cor- 
rections in Eq.(B.2) are taken into account. For inter- 
faces, the depinning transition with a nonzero Ai|| is 
thought to be equivalent to directed percolation depin- 
ning Assuming that transverse fluctuations still do 
not affect longitudinal ones, for d — 1 the critical expo- 
nents C|| and v arc related to the correlation length ex- 



ponents i^ii 



(DP) 



and 



,{DP) 



of directed percolation through 



(DP) I (DP) 



0.63, while 



1^—1^11 ~ 1.73 and (^\\ — ty ^ it^^^ 
the dynamical exponent is z\\ = 1. This in turn gives 

= - Cii )z. = ^f^)-i.p« 0.64. 

Using the connection to interface depinning further, 
we next consider tilting the FL away from the symme- 
try axis c. In this case, {fxf\\) and {fxf±) are nonzero, 
and Fc depends linearly on S||, leading to terms propor- 
tional to 9x^11 in the equation of motion. These further 
suppress the roughness exponent to C|| = 1/2 Q. The 
analysis of transverse fluctuations for these two situa- 
tions, and many other possible ones, are complicated by 
the absence of a suitable perturbative treatment. Differ- 
ent types of anisotropy may lead to distinct transverse 
exponents even while the longitudinal ones remain iden- 
tical. (Similar to the difference between Models A and 
B, although the latter is unstable to the inclusion of non- 
linear terms.) To systematically search for universality 
classes, we may start with the most general equation of 
motion, which has the gradient expansion. 



+ -^K,i3jdxri3dxr.y + fa{x, r,dxr, ■■■)-{ , 

and with force-force correlations that depend on dxJ^- De- 
pending on the presence or absence of various terms al- 
lowed by symmetries, these equations encompass many 
distinct universality classes. The cases that were dis- 
cussed so far are summarized in Table |. 



Situation 




C|l U Z f3 Cj_ Zj_ 


Anisotropic medium, 
generic direction 


Kafl 


0.5 2 1 1 ? ? 


Anisotropic, FL 
along symmetry axis 


= 

All! /O 


0.63 1.73 1 0.64 ? ? 


FL along symmetry 
axis, linear terms 
only (Model B) 


Kafi = 

<^/0 


1 1 1.3 0.3 2.3 


Isotropic medium, 
</> = (Model A) 


Kal3 = 

Aill =0 


1 1 1.3 0.3 0.5 2.3 



TABLE L Critical exponents corresponding to some of the 
universality classes associated with vector depinning. Entries 
in the first two rows are from Ref . ^l| : Transverse exponents 
are not known and these cases may correspond to more than 
one universality class identified by distinct (^±, z±. 



dtra = fJ-ai3Ff3 + Kapdxrp + Kapdlrp 



(8.4) 
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C. Generalizations 

In many systems, the dynamics involves a wide range of 
relaxation times. It is sometimes possible to average over 
"fast" degrees of freedom to obtain an effective equation 
of motion for "slow" variables. For example, the motion 
of atoms in a metal can be described by an effective the- 
ory that involves only positions of the ions, assuming that 
the electronic wavefunction always adjusts to the instan- 
taneus ionic coordinates. Similarly, the critical dynamics 
of a slowly moving solid-liquid- vapor contact line can be 
described by assuming that the liquid- vapor interface in- 
stantaneusly finds the minimum energy surface dictated 
by the position of the contact line ||3^. The elimination 
of these additional degrees of freedom may cause effec- 
tive nonlocal interactions between the remaining modes, 
which in turn acquire a different dispersion law. For ex- 
ample, in contact line dynamics, the elastic energy asso- 
ciated with a mode of wavevector q is proportional to |g| 
instead of q^. In general, one may consider a situation 
where the elastic energy is proportional to \q\'^ for some 
value of a. The scaling analysis can be easily generalized 
to such cases; the most important change is the modifi- 
cation of the upper critical dimension to dc = 2a. The 
exponents can be easily calculated for general cr, as was 
done by us for the critical dynamics of a contact line |3^] 

(^ = 1)- ... 

The possibility of experimental verification of our re- 
sults lies in the ability to accurately measure the motion 
of individual FLs and the noise spectra (for both nor- 
mal and Hall voltages) generated by FL motion. Very 



recently, there have been successful experiments that de- 
tected the thermal motion of individual FLs at nomi- 
nally zero magnetic field and bulk current using SQUID 
probes, and analyzed the noise correlation between the 
two ends of the FL js^] . A refinement of such techniques 
may eventually enable a direct comparison of theoreti- 
cal results with experiments. For example, it is known 
that the Hall angle changes sign as a function of temper- 
ature in certain superconductors |3^. It would be par- 
ticularly interesting to observe the increase in transverse 
roughness (thus the Hall Voltage noise) as the Hall angle 
approaches zero. Ultimately, it is very desirable to un- 
derstand the properties of many FLs (solid or glass) near 
depinning, especially since this situation has much more 
experimental and technological relevance. One should 
then start from a coarse-grained theory for the displace- 
ments u(x, t) of the FLs with respect to their equilibrium 
positions in the Abrikosov lattice and hope to establish 
a similar RG scheme. However, there are certainly addi- 
tional complications, such as entanglement and plas- 
ticity effects, which are difficult to incorporate in 
such an approach. 
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APPENDIX A: THE GAUSSIAN THEORY 



In this appen dix, we compute all nonzero expectation values for the Gaussian theory, described by the effective 
action 5*0 in Eq.(6.1). This is accomplished by inverting the quadratic form, as 



(R(q,c^)R^(-q,-Lj)), (R(q,tj)R^(-q,-c^)), 
(R(q,Cc;)R^(-q,-L^)), (R(q, u;)R^(-q, -co)). 



-C{uj) G-l(q,c.) 



Lct-^q, 







Gt(q,a;) 

LG(q,a;) G(q, c^)C(c^)Gt(q, c^) J 

For the case of Model A, the individual matrices are diagonal and the correlation functions can be calculated easily, 
as given in Eqs.(3.S). For the more general case of Model B, let us first consider the w ^ limit. Since R±^ occurs 
in the combination R±/v, expectation values {RaRA_) and {RaRj_) contribute at most 0{v) at the momentum-shell 
integration step. Thus, the contractions that are important for the momentum-shell integration are and 
Setting V — and inverting the matrix yields 



G{q,Lu) 
GCG^q, Lo) 



Knq^ 
K2iq^ 



ILoAll 

iu}A2i 



-iLuAl2 
-iujA22 



1 



1 



Kll'q' 



C{lo) 



(Al) 
(A2) 
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which leads to Eqs.(6.29). To determine the full form of the correlation functions in a renormalized Gaussian theory, 
we need to perform a full matrix inversion. In the small v limit we obtain 



where 



|i?ii(q,^)P> = 



detG 



Wq^ [Kl^Cw {io) - i^22i^i2ReCx (co) + Kl^Ci_ {u)] ] , 



detG-i|2 « [K^^\* + py] v^q^ {[KnK22 - K2iKu]/K^^)' + Al, 



Similarly, 



(|i?^(q,c.)p) - idetG-112 {(-^2iq|H-^2ij^iiRe[CxH] + j^i\C±H) q' 



-ujq^{KiiA2i - K2iAii)lm{C^ {lj)] 

W {Al^C\\{oj) - A2iAnRe[Cx(^^)] + A\^C^{u)) ] 



At the fixed point found for Model B, Eqs.( 6.32 ) are satisfied, and the correlation hmctions simplify to 

C{lo) ^ ( K^q^ 



(|i?ii(q,^)P> = 
(l^±(q,^)r) = 



,2,. ,2 II \ {ujIv)p\ 



where 



T\\ [x) 



K11K22 — K21K12 



1 + a; 



Kii - KK21 

f jK22^Ki2Miq_ 

\2{KiiK22 - K21K12) 



P-L 



Kn + KK21 



J_ / Kii + KK21 All - kA2i ^ 
\Kii- KK21A11 + KA21, 



Kii - KK21 

L J L cons 

r cons 
11, 



22, 



x < 1 
const, a; 3> 1. 







1 " 


)] 


1 


1 + - 









const, a; ^ 1, 
X > 1. 



(A3) 
(A4) 



The functions describe crossovers of the overall amplitudes of the correlations, due to the coupling between 
longitudinal and transverse modes. 



APPENDIX B: VERTEX RENORMALIZATION 



In this appendix, we derive recursion relations for the renormalized vertex functions = cj^\u). Let us 

start by considering Ua,n{u) for a given n. As usual, we split the fields R = and R = R^ +R^, where fields 

with the superscript ">" correspond to fluctuations within the momentum shell Ae~^^ < g < A, which are averaged 
over. In evaluating {e^)^ , we encounter two types of nonzero contractions. 



i?>(q,t)i?>(-q,r 



i?>(q,i)i?>(-q,<'; 



All 
1 



exp 



Kq^{t' ~t) 



All 



e(t' - 1) 



1 

X2A4 



Ua,Q {v(t-t')). 



within the momentum shell Ae^^^ < q < A, and for time scales t — t' ^ 0{l/v). (From now on, we suppress the 
subscript for notational simplicity.) Contributions to the renormalization of Ua,n come from both (iY)^ and {U'^)'^ , 
as 
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= E 2\in + 2)\ I ^'^'^^i^^^ C/a.„+2(l - 2) (i?„(l)i?„(2)[i?||(l) - (2)]"+2)^ + • • • 

a 

= E 2!(nV2)! C2 / '^'^'^^i^^ f/a,„+2(l - 2)i?<(l)i?<(2) 

X [i?<(l) - i?<(2r ([i?>(l) - i?>(2)p) + . . . , 
with obvious abbreviations for the arguments of U, R, R. Evaluating the expectation values, we get 

([i?>(l)-i?>(2)p) ^ 1^ {(i?>(q,ti)i?>(-q,ti)) + (i?>(q,i2)i?^(-q,t2)) 

^2^R>iq,h)R>{-q,h))} 

where denotes integration over the momentum shell and Sd is the surface area of a unit sphere in d-dimensions. 
Thus, the correction to U^j^{u) from (U)-' is equal to 

^^ifdC/a,„+2H[C/||,o(0) - t/||,o(")], 

where Kd = A'^~'^Sd/[{2TTYK'^]. The contributions from {U'^)^ are similarly calculated, as 

n+i If r 

{U')>^ — / d''xdhdt2 / dV<d4(7,,„,(l-2)C/,.„+2_„,(l'-2') 

X (i?„(l)i?„(2)i?^(l')i?7(2')[i?||(l) - i?^||(2)]"[i?^||(l') - i?||(2')]"+'-")^ + • •• 

n+l f f 

= E 2!(m-l)!2!(n+l-m)! J ^'-^hdh J " 2)C/||,„+2-™(l' - 2') 

X [R<{1) - i?<(2)]"-i[i?<(l') - i?<(2')]"+'"" 
x(i?||(l)i?||(2)i?||(lOi?||(20[4(l)-i?f(2)][i?>(lV4(^^^^ 

+2 y y , , ^ + 2 - "A /■ /■ , 

2!to!2! n + 2-m ! V 2 /i J ^ ^ 

xC/„,™(l - 2)C/i|,„+2_™(l' - 2')i?<(l')i?a (2')[i?f (1) - i?f (2)]" 

x[i?<(l') - i?<(2')]"""' (2)[i?>(l') - i?^(2')]'\ + • • • . (B2) 



The eval uati ons of the expectation values are tedious but straightforward. As an example, let us evaluate the second 
half of Eq.(p2|) explicitly. First of all, 

i?>(l)i?>(2)[i?>(l')-i?f(2')]^ 

i?>(l)i?>(l')\ (i?^(2)i?>(l')) + (i??(l)i?^(2')) (i?^(2)i?>(2') 



II v-^;^^|| ; / \-''-|| 

-2 (i?>(l)i?>(l')) (i?>(2)i?>(2')) - 2 (i?>(l)i?>(2')) (i?>(2)i?>(2') 

The first two terms do not contribute to U^J(u), since they are proportional to b{t\ — t[)6{t2 — t'l) and 6{ti — 
^2)^(^2 — ^2) respectively. (These delta functions force ti to be equal to t2- Since the expectation value is multiplied 
by — i?jj~(2)]'", the final contribution is zero.) The last two terms are equal to 

" d^q } dW , ... ^ exp{i(q + q') • (x - x')} 



w){Ke) 
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Integrating over t[, t'2, x' (which yields 5''(q + q')) and subsequently over q', the second half of Eq.(B2) becomes 



X {C/„,™(1 - 2)C/||,„+2_™(1 - 2) + (-1)"-"C/„,™(1 - 2)C/||,„+2_™(2 - 1)} . 
The first half of Eq.(B2) can be evaluated similarly. The full contribution to C^q,„(u) from {U^)^ is thus equal to 

I. m=l ^ 

+ E Q ^ [t^o,m(")C/||,„+2-™(") + {~ir+^-"'U^,MUln+2-rn{~n)]\ . 
7n—l ^ ^ 

(In the expansion of {e^)cj there is a factor of 1/2 in front of {U'^)^ .) Adding all contributions, the effective vertex 
fmiction Uf{u) is found to be 



n+l 



n + l 



(B3) 



provided that 



Uc^Mu) = (-l)'"?7„,„(-u). 



(B4) 



Under the scale transformation (6.5), which brings the momentum cutoff to its original value, we see that u 
(1 + C\\Si)u. Thus, the renormalized vertex function is given by 



Uln{u) = Uln{u)+Si 



dl 



= ^tn ((1 + C\\St)u) {1 + 5^[d + 2z|| + 2(6*11 -d) + nC||]}. 



(B5) 



Keeping only terms linear in 6t, and identifying [/|| „(m) with the nth derivative of C||(u), we finally obtain the 
differential recursion relation for C||(u): 



= [e + 20|| + 2(z|| - 2)]C\\{u) + (^^(w) - ifd {[^(u)]^ + [C\\{u) - q, (0)]q(u)} 



(B6) 



Note that the identification of [/|| „(w) with the nth derivative of C\\ {u) is sel f-co nsistent, since recursion relations for 
?7|| „(m) are correctly recovered by taking n derivatives of Eq.(B6). Also, Eq.(B4) is automatically satisfied when this 
identification is made since C\\{u) = C\\{—u). 

A similar computation can be performed for C±_ {u) , yielding 



Upon rcscaling, the renormalized vertex function is 



U^,n{u) = U^,n{u)+5i- 



d£ 



Utn ((1 + C\\Si)u) {1 + Si[d + 2z|| + 2(61^ ~d) + nC||]}. 



Thus, we obtain the recursion relation 
dC±{u) 



d£ 



[e + 26^ + 2(z|| - 2)]Cx(u) + C\\uC'Au) - i^rf {[q|(w) - C\\iO)]C'iiu)} • 



(B7) 



(B8) 



(B9) 
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APPENDIX C: HIGHER-ORDER DIAGRAMS 



In this appendix, we show that the sum of aU contributions to the renormahzation of C±{u) from the momentum 
shell integration step vanish in the limit u — > 0^. This was already explicitly demonstrated for the leading order 
contributions that come from {iU)c and (W^)c- Since the only nonzero contractions involve R\\ and all contributions 
to the renormahzation of C±{v{t — t')) due to {l^™")^ arise from terms of the form 



(W'")c =E / d'^dtdt' ^^-"^f^, ''^^ R<{t)R<{t') 



2!n! 



Y[ d'^Xi dU dt: 



2!n,! 



'\]ni 



The expectation value clearly goes to zero as {t — t')" 
in the t — > i'+ limit. This gives us the desired result 
that C_l(0) is unrenormalized to all orders in perturba- 
tion theory. 



appears in the RG contribution is identically zero due to 
the causality of the propagator. Therefore, the singular 
part of C|| (0) does not enter the renormahzation of An 
(or yO|| in the case of Model B) to one-loop order. 



APPENDIX D: HIGH-FREQUENCY 
STRUCTURE OF U« 2 



APPENDIX E: RENORMALIZATION OF MODEL 



In this appendix, we shall demonstrate that there are 
no divergences in the renormahzation of An, at least 
to 0(e). In order to do this, we examine the full form of 
the bare vertex function f7|| 2 obtained from MF theory. 



C^||,2(ii,i2;^i,i2) — 



{r\\{ti)r\\{t2)) 



MF,c 



The low- frequency analysis of this vertex function gives a 
result proportional to l/v when all times are within 0(1) 
of each other. This may potentially give an 0(1/ v) con- 
tribution to the renormahzation of An. Indeed, an ex- 
ternal impulse of magnitude e right before a "jump" (the 
fast motion between consequent local minima) shifts the 
jump time by e/77 and creates a response of 0{l/v) right 
after the jump takes place. However, an impulse right 
after a jump does not affect the jump time and creates a 
response of only 0(1). Thus a singular response is seen 
if all times are in the vicinity of a jump, say, at time tj. 
t^ll , 2(^17 ^2; ^'17*2) can be as large as (^(w"^) if t[ and t'2 are 
both slightly less than tj, and ti and ^2 are both slightly 
greater than tj . Considering that the probability of be- 
ing close to a jump is v, this term can potentially con- 
tribute as much as 0(v^^) to the renormahzation of An 
upon statistical time-averaging. A careful analysis, and 
explicit evaluation of this vertex in the case of a periodic 
potential |^^, show that this is the only way a singular- 
ity may occur in the RG contributions. However, when 
the times ^'1,^2 of fields i?|| are smaller than the times 

ti,t2 of fields the contraction (^R\\{ti)R\\(t'j)^ which 



Details of the RG calculation for Model B are pre- 
sented in this appendix. For the sake of brevity, we shall 
only consider the renormahzation of the parameters in 
the Gaussian theory, i.e. the propagator, and the two- 
point correlation functions U\\^o{u), U±,o{u), [/x.o(u). 
The renormahzation of higher-order vertex functions are 
again related to derivatives of Ca through C/a,„(u) = 

Nonzero co ntrac tions involved in the calculation are 
given in Eqs.( |6.29 ). The parameters ^4x2,^22 (thus k), 



and Ka-y (thus K\\, K^, and pi_) do not get contributions 
from the momentum shell integration, and give rise to 
exponent identities discussed in the text. On the other 
hand. An and A21 (thus as well as the functions 
Ca{u), are renormalized. Let us start by looking at the 
renormahzation of two-point correlation functions Ca{u). 
By definition, C\\{u) = C||(— u) and C_l(u) — Cj_(— u), 
but Cx(m) 7^ Cx(— m) in general. It is convenient to 
write Cx(u) in terms of its even and odd parts Cxs{u) 
and Cxa(u) respectively, and follow their renormahza- 
tion separately. 

The momentum shell integration procedure is similar 
to the one presented in Appendix albeit more cum- 
bersome due to many more nonzero contractions. Nev- 
ertheless, carrying out the computation yields 



(El) 



for 7i > 0, where 



X||(u) = q'(w)[C(w) - (7(0+)] +C[|(u)C'(w) 

{cl{u)C'Au) [C'xs(«)/2]2 + [C'x^(«)/2]2} 
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+kC^a{u) {c[,(0+) + «[C',^(0+)/2]} , 
J^iu) = C'l{u)[C(u) - C(0+)] + C'^{u)C'(u) 

- {c[{u)C'^{u) - [C'^s{n)l2Y + [C'^Au)/'^?} 
+C^a{u) {[C'^AQ+)I2]- kC'AQ+)} , 
I>^s{u) = C";s{u)[C{u) - C(0+)] + C"^si^^)C"iu) 

+2k [c[ {u)C'Au) - [C'^s{^)l2f + [C",^(u)/2]2} 

-Cxa(m) [c[ (0+) + nC'^A{Q+) - A«2ci(0+)} , 
1^a{u) = C':,a{u)[C{u) - (7(0+)] + C'^Aiu)C'{u) 
+4«{q (u)Ci(0+) -q(0+)Ci(u)} 

-C'^si^) {C'l (0+) + «C^xaH - «'C'i(0+)} . 

Thus, the renormaUzation of C{u) is given by 

C<{u) = C{u) - SlKd [c"{u)[C{u) - (7(0+)] + [C'{u)f] , 

which leads to the functional recursion relation 
dC{u) 



d£ 



+ 20|| + 2(z|| - 2)]C{u) + C||uC'(m) - Kd{[C'{u)]^ + [C{u) - C(0)](7"(u)} . 



(E2) 



(E3) 



(E4) 



(E5) 



(E6) 



(E7) 



This is identical to Eq.(p36|), with the substitution C||(u) — > C{ u). I t is straightforward to verify that there exists a 
fixed point where individual matrix elements Cq(u) satisfy Eg. ( 6.3^ ). (Cxa(u) = at this fixed point.) 

Let us next examine the renormalization of pn . Leading order contributions come from {U)q , and a calculation 
along the lines presented in Sec. VI gives 



A< = An 



(27r)'ip|| 



(27r)'^Pll 



which can be combined to yield 







(E8) 



The fixed-point function C* (w) is identical to that of C*^ (u) in Model A, and its behavior near u = is also given by 
Eq.(|^. Thus, we obtain 



p< = Pll - 6ep\\KdC"iO+) = - 6i{2e/9)] 



which leads to the recursion relation 



^=P||[0,l+C||-26/9]. 



(E9) 
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